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Abstract 

Exact analytic expressions for planetary orbits and light trajectories 
in the Schwarzschild geometry arc presented. A new parameter space 
is used to characterize all possible planetary orbits. Different regions in 
this parameter space can be associated with different characteristics of the 
orbits. The boundaries for those regions arc clearly defined. Observational 
data can be directly associated with points in the regions. A possible 
extension of these considerations with an additional parameter for the 
case of Kerr geometry is briefly discussed. 

PACS numbers: 04.20.Jb, 02.90.+p 



1 Introduction 

Nearly a century after Einstein's theory of general relativity was found to cor- 
rectly predict the precession of the planet Mercury around the sun and the 
deflection of light by the sun's gravitational field, the problem of understand- 
ing orbital trajectories around very massive objects still retains interest as it 
relates to current astrophysical topics [1] such as the study of gravitational 
waves. Among the numerous works on this subject, we should mention the 
classic publications of Whittaker [2] , Hagihara [3] and Chandrasekhar [4] on the 
Schwarzschild geometry, and the more recent work of Levin and Perez-Giz [5] 
on the Schwarzschild and Kerr geometries. 

In the work of Chandrasekhar and that of Hagihara, the orbits are classified 
into various types according to the roots of a certain cubic erjuation, while in the 
work of Levin and Perez-Giz, the orbits are classified topologically by a triplet of 
numbers that indicate the numbers of zooms, whirls and vertices. In the work 
of Levin and Perez-Giz, the orbits were obtained by numerically integrating 
the integrable equations. These authors used the planet's energy and angular 
momentum as the principal physical parameters, and made extensive use of an 



effective potential for describing the Schwarzschild orbits, as most studies on 
the topics of general relativity do. 

In this paper, we first present, in Section 2, three explicit analytic expressions 
for the orbits in the Schwarzschild geometry: one is for periodic orbits [6] (in- 
cluding a special case that we call asymptotic), and two axe terminating orbits. 
The explicit analytic expressions that we derive not only describe the precise 
features of the orbits (periodic, processing, non-periodic, terminating, etc.) but 
also clearly indicate two physical parameters which can be used to characterize 
these orbits. These two dimensionless parameters are specific combinations of 
the following physical quantities: the total energy and angular momentum of 
the planet, the masses of the massive object and the planet, and, of course, the 
universal gravitation constant G and the speed of light c. These two physical 
dimensionless parameters were first used by one of us in rcf.6. We shall refer to 
these two quantities as the energy eccentricity parameter e and the gravitational 
field parameter s respectively (or simply as the energy parameter and the field 
parameter). They will be defined in Section 2. We will use neither the common 
convention of setting G = c = 1, nor the energy and angular momentum of the 
planet by themselves, as the physical parameters for characterizing the orbits. 
With the energy parameter (0 < e < 1) plotted on the horizontal axis and the 
field parameter (0 < s < oo) plotted on the vertical axis, the parameter space 
for all possible orbits will be shown to be divisible into three sectors, which we 
call Regions I, II and IF, that have clearly defined boundaries. Region I permits 
periodic and terminating orbits. Regions II and IF terminating orbits only. 

In Section 3, we describe Region I for the orbits in greater detail. We first 
divide Region I by lines each of which represents orbits described by elliptic 
functions of the same modulus k. We then give a more physical division of 
Region I which consists of nearly horizontal lines each of which represents orbits 
that have the same precession angle A^, and of bent vertical lines each of which 
represents orbits that have the same "true" eccentricity e. The terminating 
orbits will be characterized by two parameters one of which is the angle at 
which the planet enters the center of the blackhole, and the other being the 
initial distance of the planet from the star or blackhole. In Sections 4 and 5, 
we describe Regions II and IF in which all orbits are terminating, and we again 
divide Region II by curves of constant modulus k each of which describes orbits 
with the same modulus. Regions II and IF are separated by the Schwarzschild 
horizon. Thus the grid of our "map" can be used to describe all possible orbits in 
the Schwarzschild geometry in their entirety. The observational data related to 
a planet's orbit about some giant star or blackhole can be directly identified with 
a point having certain coordinates (e, s) on our map, which can then be used for 
estimating the physical characteristics associated with the star or blackhole and 
that of the planet itself, assuming that the star or blackhole is not spinning very 
fast. For the Kerr geometry, another dimensionless quantity, which is clearly 
the ratio of the spin angular momentum per unit mass of the blackhole to the 
orbital angular momentum per unit mass of the planet, should enter into the 
consideration. In Section 6, we briefly discuss a possible extension of our results 
to the case involving a slowly spinning blackhole, at least to the first order 
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perturbation, by rcscaling the physical parameters involved. 

In Section 7, we study the deflection of light by the gravitational field of a 
very massive object. A single dimensionless parameter will be used to charac- 
terize the region. We show that here too, wc should divide the region into three 
sectors, which we again call Regions I, II and IF, and we present three analytic 
expressions for the trajectories of light applicable in these different regions. Re- 
gion I has trajectories of light that get deflected, and Regions II and IF have 
trajectories of light that are absorbed by and terminate at the blackhole. 

In Section 8, we give a summary of our results. Proofs of many interesting 
analytic relations among the parameters appearing in these studies arc given in 
several appendices. Since our results presented in this paper cover gravitational 
fields of all ranges, from the weak field produced by the sun of our solar sys- 
tem, for example, to the very strong flcld produced by a blackhole, wc want to 
avoid referring to the massive object that produces the gravitational field as a 
blackhole, and prefer to refer to it as the star or blackhole, and we shall refer 
to the object of a much smaller mass that orbits around it as the planet or the 
particle. 

We have supplemented our many analytic results with numerous tables that 
present various physical quantities such as the minimum and maximum distances 
of the planet from the star and the angles of precession of the orbits that are 
calculated from our analytic expressions, as well as numerous figures that show 
various kinds of orbits of the planet and various kinds of defiection of light. 



2 Analytic Expressions for the Orbits 

We consider the Schwarzschild geometry, i.e. the static spherically symmetric 
gravitational field in the empty space surrounding some massive spherical object 
such as a star or a blackhole of mass M. The Schwarzschild metric for the empty 
spacetime outside a spherical body in the spherical coordinates r,6,^ is [1] 

df = c2 (l - ^) df - (l - ^) dr^ - r^de^ - sin^ edcj>^ (1) 



where 



2GM 

<^=—^ (2) 



is known as the Schwarzschild radius, G is the universal gravitation constant, 
and c is the speed of light. If [x'^] = {t, r, 6, (p), then the worldline x^(t), where 
T is the proper time along the path, of a particle moving in the equatorial plane 
6 = 7r/2, satisfies the equations [1] 



{l-^)t=n, (3) 
-(l-^y\^ -r^cP =c\ (4) 
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r^cj> = h, (5) 

where the derivative • represents d/dr. The constant h is identified as the 
angular momentum per unit mass of the planet, and the constant k is identified 
to be 

E 

K = 



1 ' 



where E is the total energy of the planet in its orbit and toq is the rest 
mass of the planet at r = oo. Substituting eqs.(3) and (5) into (4) gives the 
'combined' energy equation [1] 



r 

Substituting dr/dr = {dr / d(f){d(j)/ dr) = {h/r'^){dr /d(j>) into the combined 



energy equation gives the differential equation for the orbit of the planet 

au^ -u^ + Bu + C (7) 



du\'^ 

d4>) 

where u = 1/r, B = 2GM/h^, C = c^{k^ - l)/h^. Following Whittaker [2], 
it is convenient to change variable from u to a dimensionless quantity U defined 

by 

u=l(^~D^Uau-D. (8) 



or u = 4:11 1 a + l/(3a) so that eq.(7) becomes 

= 4C/3 - g^V - 33 (9) 



where 



Qi = s 

12 

03 = s H — (1 — e IS , 

216 12 4^ ' ' 



(10) 



and where 



and 



h^c^jn^ - 1) 
(GAf)2 



1/2 



(11) 



GM 
he 



(12) 
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The two dinicnsionlcss parameters e and s, which are defined by the two 
above equations and which we call the energy and field parameters respectively, 
will be the principal parameters we shall use for characterizing the orbit of a 
planet. It will be noted that the constant c^{k'^ — 1) which is < for a bound 
orbit, can be identified with 2EQ/m in the Newtonian limit, where Eq is the 
sum of the kinetic and potential energies and is < for a bound orbit, and m 
is the mass of the planet (which approaches mo), and that 



m{GMY_ 



(13) 



is the "eccentricity" of the orbit. Also, in the small s limit, the orbit equation 
can be shown to be given by 

^^^[l-ecos(l-5)0], (14) 

where 5 ~ 3{GM)'^/{hc)'^. Thus r assumes the same value when (p increases 

to (f> + 27r/(l — 5). Comparing this with the increase of (f> from (/) to (j) + Iti. the 
ellipse will rotate about the focus by an amount which is the angle of precession 

27r „ „ , 67r(GM)2 

- - 27r ~ 1-Kb = . 15 

1 — h^c^ 

This is the well known approximate expression for the precession angle for 
the case of very small s. The limiting case for 5 = is the well known orbit 
equation in Newtonian mechanics. We should note that while the limit s = 
(and thus ^ = 0) cannot be strictly correct in principle so long as M 7^ 0, this 
limit can be used for many practical cases with great accuracy as evidenced 
by the predictions of Newtonian mechanics. A special case of these Newtonian 
orbits is the circular orbit of radius r = h? /GM for e = 0. 

We now derive the exact analytic solutions of eq.(9) and classify the three 
possible solutions from a purely mathematical viewpoint, and we shall consider 
their physical interpretations in the next section. We first define the discrimi- 
nant A of the cubic equation 

4[73 _g^u-gs = Q (16) 

by 

A = 27.g2-.g3. (17) 

The three roots of the cubic equation (16) are all real for the case A < 0. 
We call the three roots 61,62,63 and arrange them so that 61 > 62 > 63; the 
special cases when two of the roots are equal will be considered also. For the 
case A > 0, the cubic equation (16) has one real root and two roots that are 
complex conjugates. The analytic solutions of eq.(9) that we shall present will 
give the distance r of the planet from the star or blackhole in terms of the 
Jacobian elliptic functions that have the polar angle (j) in their argument and 



5 



that arc associated with a modulus k that will be defined. Instead of writing 
r, we use the dimensionless distance q measured in units of the Schwarzschild 
radius a and defined by 

r 1 

q^-^—. (18) 
a au 

The dimensionless distance q is related to U of eq.(8) by 

- = \+4U. (19) 
q 3 

We now give the three analytic solutions of eq. (9) in the following. 
Solution (Al) For A < 0, ei > 62 > C/ > 63. 

Writing the right-hand side of cq.(9) as 4(ei — U){e2 — U){U — 63), cq.(9) 
can be integrated with <p expressed in terms of the inverse Jacobian sn function 
[7]. After a little algebra and some rearrangement, the equation for the orbit is 
found to be 



- = ^ + 463 + 4(62- 63)572^(7^, fc) 

q o 

= 3 + 463 + 4(6.-63) ^^^^^;;;^; , (20) 

where the point at (f> = has been chosen to give U = e^. The constant 7 
appearing in the argument, and the modulus fc, of the Jacobian elliptic functions 
are given in terms of the three roots of the cubic equation (16) by 



7 = (el-e3)^/^ (21) 

^2 ^ (22) 

61 - 63 



where 61, 62, 63 are given by 
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Equation (20) was first given in rcf.6 using a slightly different approach that 
was initiated by Whittaker [2]. In addition, eq.(20) was shown to reduce to 
eq.(14) for the case of very small s which in turn gave the known approximate 
precession angle given by cq.(15). The modulus k of the elliptic functions has 
a range < < 1. Since the elliptic functions sn, cn and dn are all periodic 
functions of ^ for < k"^ < 1, we shall call this solution for the orbit the 
periodic solution. For the special case of fc^ = 1, since sn{j(f>,l) = ta.nh{'y(f>) , 
cn(7(/>, 1) = dn{'~^4',l) = sec/i(7^), the solution is no longer periodic, and we 
shall refer to it as the asymptotic periodic solution. 

The period of cn{2-f(f),k) is 4:K{k), and the period of dn{2^4),k) and of 
sn^(70, /c) is 2K{k), where K{k) is the complete elliptic integral of the first 
kind [7]. For fc = 0, sn{x,0) = sinx, cn{x,0) = cosx, dn{x,0) = 1. As fc^ 
increases from to 1, K{k) increases from tt/2 to oo. The distance r of the 
planet from the center of the star or blackhole assumes the same value when its 
polar angle (f) increases from (f) to <p + AK/{2'y) = (p + 2K/'y. Comparing this 
with the increase of (j) from (/> to + 27r in one revolution, the perihelion (or the 
aphelion) will rotate by an amount 

A<^=^^-27r, (25) 

7 

which is the exact expression for the precession angle. For k'^ close to the 
value 1, the planet can make many revolutions around the star or blackhole 
before assuming a distance equal to its initial distance. Thus if n is the largest 
integer for which 2K{k)/^ is equal to or greater than 2n7r, the angle of preces- 
sion should be more appropriately defined as 2K(k)/^ — 2mr. For the sake of 
consistency, however, we shall stick to the definition given by eq.(25). 

For the case of very small s and to the order of s^, it was shown in ref.6 that 
7 ~ [1 - (3 - e)s2]/2, fc2 ~ 4es2^ K{k) ~ 7r(l + es'^)/2, and substituting these 
into eq.(25) gives the well known approximate result given by eq.(15). 

For these periodic orbits, we note that the maximum distance rmax (the 
aphelion) of the planet from the star or blackhole and the minimum distance rmin 
(the perihelion) of the planet from the star or blackhole, or their corresponding 
dimensionless forms Qmax (= ^'max/a) and qmin {— ^min/ct), are obtained from 
eq.(20) when jcf) = and when 70 = K{k) respectively, and they are given by 

-^ = ^+463, (26) 

ymax *J 

and 

— = ^ + 462, (27) 

where 62 and 63 are determined from eqs.(23), (24) and (10) in terms of e 
and s. 

Although we call the orbits given by this solution for < fc^ < 1 periodic, 
they are not necessarily closed orbits. It is seen from the precession discussed 
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above that for Ai/) = /tt, unless / is a rational number, the orbit will not close 
and it is not strictly a closed periodic orbit. However, for all practical purposes, 
any irrational number when truncated becomes a rational number, and thus the 
orbit will be closed. The distinction of closed and non-closed orbits depending 
on whether / is rational or irrational is of course of profound theoretical interest 
[5]. 

For a general periodic orbit that precesses, the general or true eccentricity e 
of the orbit is defined by 

^max ^min Qmax Qmin tc\r,\ 

£= = , (28) 

^max ~t~ ^min Q'max ~t~ Qmin 

where (/max and gmin are given by eqs.(26) and (27). 

We shall show in the following section that the true eccentricity e is in general 
not equal to the energy eccentricity parameter e defined by eq.(ll), but that 
£ ^ e in the limit of s ^ 0, i.e. in the Newtonian limit. For the special case of 
£ = 1 however, we shall show that it coincides with the special case of e = 1 for 
all applicable values of s, and that it signifies an unbounded orbit. 

We now proceed to present the second solution. 
Solution (A2) For A < 0, i7 > d > 62 > 63. 

We write the right-hand side of eq.(9) as 4({7 — e\){U — e2){U — 63) and 
eq.(9) can be integrated with (f) expressed in term of the inverse Jacobian sn 
function [7]. The equation for the orbit is found to be 



4^^^^i7rf4F^' (29) 



1 ^ 1 _^ ^ ei - e2sn^(70, k) 
q 3 cn^(70, fc) 

where 7, fc, ei, 62 and 63 are given by eqs.(21)-(24) as in the first solution. 
This solution gives a terminating orbit. The point at ^ = has been chosen to 
be given by 

- = ^ + 4ei. (30) 

The planet, starting from the polar angle </> = at a distance gi from the 
blackhole, plunges into the center of the blackhole when its polar angle 4>i is 
given by cn(7^i, fc) = 0, i.e. when 

7 

where 7 and k are given by eqs.(21) and (22). 

The region of (e, s) where orbits given by solutions Al and A2 are applicable 
will be called Region I, and it will be described in greater detail in Section 3. 
Thus each point (e, s) of parameter space in Region I represents two distinct 
orbits, one periodic and one terminating. At the same coordinate point, the 
characteristic quantities that describe the two distinct orbits are related. For 
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example, by noting 61 + 62 + 63 = and from eqs.(26) and (27), qi can be 
expressed as 



Ql V^min ymax 

where (/min and ^max are the minimum and maximum distances for the pe- 
riodic orbit at the same coordinate points (e.s). It will be noted that qi is less 
than gmin, i-e. for the terminating orbit the planet is assumed initially to be 
closer to the blackhole than the gmin for the associated periodic orbit, except at 
k"^ = I where qi = g,„in and the planet has a circular instead of a terminating 
orbit that will be explained later. 

We note that the terminating orbit equation (29) presented has no singularity 
at the Schwarzschild horizon q = I, bcc;ausc, as is well known, q = 1 is a 
coordinate singularity and not a physical singularity. The orbit obtained from 
continuing (f) beyond the value (pi = K{k)/j may become interesting if the 
concept of " whiteholc" turns out to be of physical relevance. 

For now, the orbits in Region I are characterized mathematically by A < 0. 

We now present the third solution. 

Solution (B) For A > 0. 

Define 





B = o 53-\7^ , (32) 



where gs and A are defined by eqs.(lO) and (17). The real root of the cubic 
equation (16) is given by 

a = A + B (33) 

and the two complex conjugate roots b and b are —{A+B) /2±{A — B)y/3i/2. 
We further define 

7= [3(A2^AB + B2)]1/4 (34) 

and 

fc2 ^ 1 _ HA + B) ^ 1 _ ^ .35. 

2 4^3(^2 + AB + £2) 2 472' ^ ' 

Writing the right-hand side of eq.(9), with U > a, as A{U — a){U — b){U — b) , 
eq.(9) can be integrated with </> expressed in terms of the inverse Jacobian cn 
function [7]. We find the equation for the orbit to be 
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1^1 7' + a- (7'-a)cn(27^,fc) 

g 3 1 + cn(27<?i, k) ' ^ ^ 

This is a terminating orbit. The initial distance q2 of the planet at = 
has been chosen to be given by 

- = ^ + 4a. (37) 

It plunges into the center of the blackhole when its polar angle (p — (j)2 is 
given by 

K{k) 

where 7 and k are given by eqs.(34) and (35). Again, we note that the orbit 
equation (36) has no singularity at g = 1. 

The region of (e, s) where orbits given by eq.(36) are applicable will be 
divided into two sectors called Regions II and 11', the boundary between which 
will be defined later. They have terminating orbits only. For now, the orbits in 
Regions II and IF are characterized mathematically by A > 0. 

As for the initial points of the orbits discussed above, by comparing eq.(19) 
with the orbit equations (20), (29) and (36), and with eqs.(26), (30) and (37), 
we already noted that our choice of = in our orbit equations is such that 
it gives U = 63, ei, and a respectively that in turn give q = gmax, 91, and 92 
as the initial distances of the planet from the star or blackhole. We then note 
from eq.(9) that dU/d^ = and hence dr/d<p = for the planet at these initial 
points of the trajectories, i.e. the trajectory or more precisely the tangent to 
the trajectory at = is perpendicular to the line joining the planet to the 
star or blackhole. All this will be seen in the figures presented later, and all 
our references to the initial position of the planet from here onward assume 
that the trajectory (as (p increases from 0) of the planet at its initial position is 
perpendicular to the line joining the planet to the star or blackhole. 



3 Region I 

Consider the orbits expressed by eqs.(20) and (29) given by solutions Al and 

A2 and characterized mathematically by A < 0. We call the region covered by 
the associated range of values for (e, s) Region I. 

To gain a preliminary perspective, consider the earth (as the planet) and the 
sun (as the star) in our solar system. Substituting the mass of the sun M = 
Ms = 1.99 X and the angular momentum of the earth per unit mass of 

the earth h = 4.48 x lO^^m^/s, we find s = 0.983 x 10"". The energy eccentricity 
parameter e, which is equal to the true eccentricity e of the earth's orbit for 
such a very small s value, is known to be about 0.017. The approximate relation 
k"^ ~ 4es2 gives the squared modulus of the elliptic functions that describe the 
earth's orbit to be k^ = 0.657 x 10~®. We see that for the planetary system that 
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is familiar to us, the values of s and k'^ are very small indeed. We may also note 
that the Schwarzschild radius a = 2GMs/c^ — 3 km would be well inside the 
sun which has a radius of 6.4 x lO^fcm. The earth's dimensionless distance is 
g ~ 5 X 10^ from the sun's center. For this value of ,s. with g,nin ~ gmax — 5 x 10^, 
the orbit given by eq.(29) from solution A2 would require the initial position qi 
of a planet to be ~ 1 according to eq.(31), i.e. the planet would have to be at 
a distance equal to the Schwarzschild radius from the center of the sun for it 
to have a terminating orbit which plunges to the center of the sun. Therefore 
the terminating orbit given by eq. (29) is inapplicable for our solar system. The 
periodic orbits, on the other hand, are perfectly valid. 

However, for cases when the massive object is a gigantic mass concentrated 
in a small radius such as a blackhole, all the possibilities presented here may 
arise. As the field parameter s increases from 0, the modulus k of the elliptic 
functions that describe the planet's orbits also increases. From eqs.(21)-(24), it 
is seen that several steps are needed to relate fc^ to e and s. In Appendix A, we 
show that a direct relationship between fc^ and e and s can be established, and 
it is given by 



l-18s^ + 54(l-e^)s^ _ (2-fc^)(l + fc^)(l-2fc^) 

(l-12s2)3/2 - 2(1 - fc2 + fc4)3/2 i^Sj 

= cos 61. (39) 

The C0&6 of eq.(39) is the same cosO that appears in eq.(24), and, in partic- 
ular, it is equal to 1, 0, —1 for fc^ = 0, 1/2, 1 respectively. 

The curve represented by fc^ = 1, after setting cos^ = — 1 in eq.(39), can be 
readily shown to give a quadratic equation 27(1 — e^)'^s^ — 2(1 — 96^)5^ — = 
that gives 

2 _ 1 - 96^ + v/(l - 9e2)2 + 27e2(l - e^f 
~ 27(1 -e2)2 ' ^^^^ 

for < e < 1, and = I/IQ for e = 1. Equation (40) representing k'^ = 1 
gives the upper boundary (for the values of s) of Region I (the uppermost 
heavy solid line in Fig.l); it extends from si = ■\/2/27 = 0.272166 for e = to 
Si = 1/4 = 0.250000 for e = 1, i.e. a line that is nearly parallel to the e-axis. 
Thus Region I is a region bounded by < e < 1, and by < s < si where si 
is given by eq.(40), in which the squared modulus of the elliptic functions that 
describe the orbits cover the entire range < k^ < 1. 

We now use eq.(38) to give a plot of lines of constant fc^ = 0.001, 0.01, 0.1, 0.3, 
as shown in Fig.l. These lines conveniently divide Region I into regions of in- 
creasing field strengths as k^ increases from to 1. On a point representing a 
particular fc^ and a particular e value, s can be determined from eq.(38) and 
the orbit is then given by eq.(20) using eqs.(98), (10) and (21). The values of s 
on these constant fc^ lines for the values of e = 0.1, 0.2, 1.0 are given in Table 
1 which thus give the coordinates (e, s) of the points on the lines representing 
different values of k"^. These coordinate points (e, s) from Table 1 are used to 
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give the following tables: Tables 2 and 3 give the values of (/max and qmin for the 
orbits obtained from eqs.(26) and (27). Note that the dimensionless distance q 
is in units of the Schwarzschild radius a which depends on the mass M of the 
star or blackhole corresponding to that particular coordinate point, and thus 
one should not compare q at two different coordinate points just by their abso- 
lute values alone. Table 4 presents the values of the precession angle in units 
of TT, i.e. Aip/ir, obtained from eq.(25). Table 5 presents the values of the true 
eccentricity e obtained from eq.(28). Tables 2-5 are to be used in conjunction 
with Table 1 for identifying the locations (e, s) of the corresponding quantities 
that are presented. The physical quantities presented in Tables 2-5 together 
with the coordinates (e, s) given in Table 1 now give all possible periodic orbits 
in the Schwarzschild geometry in its entirety. That is, the coordinates (e, s) of 
a planet orbiting a non-spinning blackhole can be identified if the observation 
data on rminj ?'max, £ and Acj) can be collected. Region I shown in Fig.l is where 
orbits given by eqs.(20) and (29) apply. In Sections 4 and 5, we shall discuss 
Regions II and IF which are shown above Region I in Fig. 2 where orbits given 
by eq. (36) apply. As an example of application of Tables 1-5, from the second 
row and second column of Tables 1-5 and using only two significant figures, 
for orbits with e = 0.10, s = 0.11, fc^ = 0.010, we find from Tables 2-5 that 
amax = 50,5i„m = 34, A(/>/7r = 0.079 or Acj) = 14°, and e = 0.19, i.e. orbits 
with those seemingly small values of s and fc^ give a precession angle of 14° per 
revolution that is already very large compared to those encountered in our solar 
system for which the precession angle is only 3.8" per century for the earth's 
orbit (for which s ~ 0.983 x 10"^, fc^ ~ 0.657 x 10"^, e ~ e ~ 0.017), and the 
value of the true eccentricity e of these orbits is already quite different from 
their energy parameter e. We thus appreciate that the range of values for s 
given by < s < si for Region I, where si ranges from 0.276166 for e = to 
0.25 for e = 1. is not as small as it seems (noting also that < fc^ < 1), and 
that the classical Newtonian orbits are restricted to a very narrow strip of the 
region indeed for which s ~ and fc^ ~ 0, and for which e~eforO<e<l. 

Although the lines of constant fc^ in Region I conveniently associate the orbits 
with the orbit equations for the periodic and terminating orbits given by eqs.(20) 
and (29) and with the physical parameters given in Tables 2-5, the precession 
angle A(j) and the true eccentricity e are more physically meaningful parameters 
that can be associated with the description of the orbit. The expressions for 
A(j) and e in terms of k and s are given by eq.(99) in Appendix A and eq.(lOO) 
in Appendix B. For a given value of A(t> and of e, we can use eqs.(99) and (38) 
to solve for s (and k) (using a numerical program such as FSOLVE in MAPLE) 
and thus locate its coordinate (e, s); and similarly for a given value of e and of e, 
we can use eqs.(lOO) and (38) to solve for s (and k). The relationship of e and s 
with £ is simpler for k^ = 1 and will be discussed later (see eqs.(49) and (50)). In 
Fig. 3, we present lines of constant A^/tt (that are nearly horizontal) and lines 
of constant e (that are bent vertical) in Region I, and the corresponding tables 
for their coordinates are presented in Tables 6 and 7. We note that because 
A(p given by eq.(25) depends on 7 given by eq.(21) as well as on K{k), the 
line of constant Acj) does not coincide with the line of constant fc^ except for 
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= \. Wc note also that the hne of constant e does not coincide with the 
(vertical) line of constant e except for e = e = 1. We show in Appendix B that 
it is only for a very thin strip of region, where s is between zero and some very 
small positive value, that e ^ e which applies in the Newtonian limit. Wc also 
show in Appendix B that e = e when e = 1 exactly. The distinction between 
e defined by eq.(ll) or eq.(13) with e defined by eq.(28) in the Newtonian or 
non-Newtonian theory has never been clearly recognized previously. 

With Fig. 3 which has curves of constant A^/tt and constant e in place. 
Region I is now partitioned into cells with the coordinate points specified by 
(A(^/7r,e). We have a clear idea what the orbits of a planet would be like at 
points within each cell in terms of their precession angle and true eccentric- 
ity, and the coordinates of these orbits (e, s) then give the energy and field 
parameters corresponding to these orbits. In Fig. 4, wc present examples of pe- 
riodic and unbounded orbits, plotted in polar coordinates {q,4>), corresponding 
to various precession angles of7r/6,7r/3,7r/2,7r,37r/2,27r,oo (vertically from top 
to bottom) for values of e = 0, 0.5, 1 (horizontally from left to right), where the 
star or blackhole is located at the origin. We first note that the orbits for which 
e < 1 are periodic and closed because / is a rational number in Acf) = fir for 
each one of them. The precession angle can be seen from the heavy solid line 
that marks the trajectory (as (j) increases) from the initial point at (/) = to 
the first point at which the distance from the origin is equal to the distance at 
(j) = 0. The true eccentricity of the orbits is s given by eq.(28). For example, for 
the orbit of Fig. 4 (al) for A(j) = 7r/6, e = 0, e is far from zero which can be seen 
from the gmin and grmax in the figure, and it can be more accurately calculated to 
be equal to 0.22629. For each of the unbounded orbits characterized by e = 1, 
the incoming trajectory coming from infinity at (/) = makes an angle with 
the outgoing trajectory going to infinity given by </> = 2K{k)/j = 2'7r-|- A0 from 
eq.(25), as can be seen in some of the figures presented. The case Aip/ir = oo 
corresponding to the special case of A:^ = 1 will be discussed later in this section 
for which the planet starting from (/max ends up circling the blackhole with a 
radius that approaches q^i^ (sec Fig.6d). 

Generally, if we are given a coordinate point in Fig. 3, for example, a point 
on e = 0.5 just above the Ai^/tt = 1/3 line slightly to the left of the e = 0.6 
curve (where e = 0.581431... and s = 0.194229...), then we find Ac/) = 60.4706... 
degrees or A(j)/n = 0.33594..., and part of the orbit is shown in Fig. 5. Whether 
the orbit will close on itself depends on whether A^/tt is or is not a rational 
number in principle, even though, as we mentioned before, a truncated number 
in practice is always a rational number and the orbit will be a closed one. We 
only show part of the orbit in Fig. 5 as the subsequent path is clear from the 
angle of precession and true eccentricity of the orbit and we are not concerned 
with how many "leaves" the orbit is going to create. Figure 3 (or one with even 
more curves of constant Ac^/tt and constant e) is a very useful map that can be 
used fruitfully with any observation data that arc obtained for any planet. 

Besides the special case — 1, the case of k^ = 1/2 is also somewhat special 
in that it allows many relationships to be expressed simply and explicitly. We 
present some of these simple relations for fc^ = 1/2 in Appendix C. It is to be 
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noted from Fig.l that the Hne of constant fc^ = 1/2 is very close to the boundary 
given by fc^ = 1. The line of constant = 1/2 for Region II, on the other hand, 
is closer to dividing the region approximately into two halves, as shown in Fig. 2. 
The fc^ = 1/2 curve for Region II will be discussed in Section 4. 

The terminating orbits in Region I given by eq. (29) can be characterized by 
the planet's initial position qi given by eq.(31), and by the angle (p^ at which 
the planet enters the center of the blackholc. It is interesting to note that 
even for these terminating orbits, the precession angle still has an "extended" 
meaning and use that we shall describe. It is clear from eq.(29) that the orbit 
terminates, i.e. q becomes zero when = K{k), but if the orbit is continued 
(by continuing to increase </>), q would assume its initial value at = when 
7(^' = 2K{k), producing a "precession angle" of = (j) — 2-k = 2K{k)/') — 2it 
which is equal to the precession angle for the corresponding periodic orbit at 
the same coordinate point (e, s). Since </>' = 2(pi, the polar angle at which 
the path of the terminating orbit enters the center of the blackhole is related 
simply to the precession angle of the periodic orbit by (pi = A^/2 + w, or 

TT 2 TT 

As (t>i/TT can be easily calculated from Acp/n for the periodic orbits given 
in Table 4, we do not tabulate it separately. The values of qi are presented in 
Table 8, and we note the small range 1 < (?i < 2.25 for the entire Region I. 
Examples of these terminating orbits are presented in Fig. 6. The dotted line 
represents the continuation of the orbit when (f) is continued beyond (p^. 

Before we discuss Regions II and IF. we want to describe three special cases: 
the case of fc-^ = which, as we shall see, is not of any interest but must be 
included for completeness; the case of fc^ = 1 which gives the upper boundary of 
Region I (and lower boundary of Region II) ; and the case of e = 1 which gives 
the right boundary of Region I (and of Regions II and IF) (see Figs.l and 2). 

(i) The Special Case of fc^ = 

The line of A:^ = coincides with the ,s = axis in Fig.l. To show this, we 
note that fc^ = implies 9 = from eq.(95). Substituting 9 = into eq.(24) 
gives s = when we use the expressions in eq.(lO) for §2 and 53. The we find 
g2 = 1/12 and 53 = 1/216, and from eq.(98), we find 

1 

= 6' 

62 = 63 = -^. 

Equation (20) then gives l/g = Oorg = oo, i.e. it is the limiting case of 
zero gravitational field. As we pointed out earlier, the classical Newtonian case 
is given by only a very narrow strip represented by fc^ ~ and s ~ for which 
q is large but finite. 

(ii) The Special Case of fc^ = 1 
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It follows from cqs.(38) and (39) that on the Une of fc^ = 1, cos^ = — 1. 
Thus from eqs.(24) and (17), we have 



A = (41) 

which can be identified as the "boundary" between Solutions A and B in 
Section 2. The range of s values for A = is 0.25 <s< 0.272166 for 1 > e > 
(see the discussion below eq.(40)), and for that range of s values, s < l/2\/3 = 
0.288675 or < 1/12 and therefore g2 > (see eq.(lO)). From eq.(41), the 
relation between 52 and 53 can be more precisely expressed as 




after noting that 53 is negative and 52 is positive for the values of s along 
the line fc^ = 1. Also from eq.(98), we note that 



/.92 

= "^ = Vl2' 
es = -M (42) 



The periodic orbits given by eq.(20) become 



1 = 1+2 /g2' l-5secfe(27(/)) 

q 3 V 3 l + sec/i(27^) ' ' 



Q \ 1/4 



where 



and where the values of 52 (and 33) are those given by the values of e and s 
on the line = I that are obtained from eq.(40). The orbit is not a periodic 
orbit; it is what we call an asymptotic periodic orbit. The planet starts from 
an initial position ^max at = given by 

1 l+4e3=^-4,/f (45) 



^max 3 3 

and ends up at (/) = 00 circling the star or blackhole with a radius that 
asymptotically approaches gmin given by 

(46) 

Qmin ^ 6 \ 6 

Equations (43)- (46) are explicit and simple equations that give the orbit 
equation, ^max, and g'min for k"^ = 1. In particular, it is seen from Table 3, for 
example, that g'min ranges from 2 for e = to 9/4 = 2.25 for e = 1, i.e. qmin is 
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still no less than twice the Schwarzschild radius for the strongest gravitational 
field that permits the periodic orbits. However, it is a very small number indeed 
compared to, say, Qmin — 5 x 10^ for the earth's orbit around the sun. 

On this upper boundary fc^ = 1 of Region I, the terminating orbit given by 
eq.(29) from Solution A2 becomes a circular orbit with a radius Qc = qi, where 
qi is the initial distance of the planet from the star or blackhole given by eq.(30). 
From eqs.(30) and (31) and noting that ei = 62 for fc^ = 1, we find that 

qc = qi = 5miii (47) 

given by eq.(46) (see Tables 3 and 8 for fc^ = 1). We shall refer to the orbits 

given by eqs.(43) and (47) as the asymptotic periodic and the asymptotic termi- 
nating orbits respectively of Region I. Thus the special cases given by eqs.(43) 
and (47) for fc^ = 1 of the periodic and terminating orbits given by eqs.(20) 
and (29) for solutions Al and A2 respectively clearly exhibit completely differ- 
ent behaviors from their counterparts for < fc^ < 1. Examples of asymptotic 
periodic orbits are shown in Fig.4g. Asymptotic terminating orbits are simply 
circles of radius equal to qi, as shown in Fig. 6(d). 

Using eqs.(28), (42), (45) and (46), for = 1 the true eccentricity e can be 
shown to be expressible in terms of 52 by 

.= ^^, (48) 
1-3^/52/3 

which can be solved to give s in terms of s, and then e in terms of s using 
eq.(38). We find that the coordinates (e, s) of a given 0.6 < £ < 1 on the line 
= 1 are given by 

_ V(l + s)(-3 + 5.) 

(3-£) ' ^ ' 

and 

2(3 -l-e) ^ ' 

It is interesting that eqs.(49) and (50) can be used in place of eq.(40) as 
parametric; equations for determining the coordinates (e. ,s) of the line k^ = 1 
as £ takes the values from 0.6 to 1. In particular, eqs.(49) and (50) allow us 
to see that the e = const, curves are not vertical (except for £ = e = 1), and 
they intersect the upper boundary si of Region I for 0.6 < £ < 1 (see Fig. 3). 
The £ = 0.6 curve, the boundary curve Si, and the s-axis axe concurrent at 
e = 0, s = ^2/27. We can conclude that periodic orbits with e = have true 
eccentricity in the range < £ < 0.6. As another example it can be shown 
using eqs.(49) and (50) that periodic orbits with e < 3v^/ll = 0.609836 have 
£ < 0.8. 

(iii) The Special Case of e = 1 
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The maximum or boundary value for e which is e = 1 is also a special case of 
interest. From its definition given by eq.(ll), since < 1, e cannot be greater 
than 1. In Appendix B, we show that e = 1 always gives an unbounded orbit for 
the periodic and the asymptotic periodic orbits of Region I and the terminating 
orbits of Region II, but not the asymptotic terminating orbit of Region I which 
has a radius given by eq.(47) independent of e and thus is not an unbounded 
orbit. Thus e = 1 is the boundary for e for Region I (as well as for Region II). 
Many explicitly simple relationships among s, /c, gmin, 9i; etc. have been found 
on the boundary line e = 1, and they are given and proved in Appendix B. In 
particular, we have, on e = 1 in Region I, that 

= 4(1 + F)2' 



I \ 1/2 



4(1 + A;2) 
l + fc2 



fc2 ' 



(52) 
(53) 



and 



qi = l + k\ (54) 

Examples of unbounded orbits for e = 1 are shown in Fig.4 (a3-g3). 
We shall now describe Regions II and IF for the orbit equation (36) given 
by solution B for the case A > 0. 



4 Region II 

Consider the orbits expressed by cq.(36) given by solution B and c;haractcrizcd 
mathematically by A > 0. The associated values for (e, s) in this case satisfy 
s > Si , where Si is the upper boundary of Region I given by eq. (40) . This region 
of parameter space defined by s > ,si can be naturally divided into two sectors 
which we call Region II and IF with Region II bordering Region I (see Fig. 2). 
The boundary between Regions II and IF is determined by the Schwarzschild 
radius in a manner to be described later in this section. 

We first want to prove that the lower boundary (for s) of Region II, charac- 
terized by A = as it is for the upper boundary of Region I, also gives fc^ = 1, 
where fc^ is calculated from eq.(35) for Solution B (In Section 3. we showed 
that for fc^ calculated from eq.(22) for Solution A, A;^ = 1 implies A = 0). 
Substituting A = into eq.(32) gives 



A = B = 




After noting that A{= B) is a negative value for the range of s values for 
A = 0, substituting the above into eq.(35) gives fc^ = 1. We also find from 
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eq.(34) that 7^ = = •\/3^/2 which agrees with the 7 given by eq.(44), 

and we find from eq.(33) that 



Substituting these into cq.(36) gives the same orbit equation (43) for the 
terminating orbit in Region II on its lower boundary as that for the asymptotic 
periodic orbit in Region I on its upper boundary. Thus on the boundary fc^ = 1 
the equation for the orbits in Region II does not represent a terminating orbit 
but is the same as the asymptotic periodic orbit for Region I given by cq.(43) 
(see Fig. 4, gl-g3). Also, from eqs.(42) and (55), we see that the smallest root 
in eq.(23) in Solution (A) is identified with the real root given by eq.(33) of 
Solution (B), i.e. 63 = a. Thus from eqs.(26) and (37), (72 = 9max when k"^ = 1, 
i.e. the initial distance 92 of the terminating orbit in Region II can be identified 
as the continuation of 5max of the periodic orbit from Region I. On the boundary 
of Regions I and II, the two other real roots ei = 62 given by eq.(42) of the 
cubic equation (16) agree with b = b given below eq.(33). The line — 1 
defined by eq.(40) is the boundary between Regions I and II; it is the upper 
boundary for Region I and is the lower boundary for Region II (see Fig.2). The 
above discussion also illustrates the transition that takes place: from a periodic 
orbit to an asymptotic periodic orbit to a terminating orbit, as one crosses the 
boundary from Region I to II. 

We now consider the upper boundary of Region II. We define this boundary 
to be that obtained by requiring the planet's initial position to be just at the 
Schwarzschild horizon, i.e. that oljtained by setting q = 1 initially at iji = 0. 
Setting g = 1 in eq.(36) for ^ = which is l/<z = 1/3 + 4a, we require a = 1/6, 
where a is the real root of the cubic equation (16). We then use the equation 



and substitute the expressions for (72 and 53 given in eq.(lO) into eq.(56) and 
solve for s. We find 



which we shall use as the equation for the upper boundary of Region II, for 
< e < 1. Thus Region II is a region bounded between e = and e = 1, and 
between s\ given by eq.(40) (the lower heavy solid line in Fig.2) and S2 given by 
eq.(57) (the upper heavy solid line in Fig.2), i.e. si < s < 82- The region defined 
by S2 < s < 00 and bounded between e = and e = 1 will be called Region IF, 
for which the planet's initial position ranges from just inside the Schwarzschild 
horizon up to the center of the blackhole. Since the same terminating orbit 
equation (36) applies in Regions II and IF, the division into two regions may 
seem unnecessary. However, the Schwarzschild radius is of physical significance, 
and it is useful to know the location of the curve S2 in the (e, s) plot which 




(55) 




(56) 
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indicates that the initial position of the planet is at the Schwarzschild horizon. 
Separating out Region IF also makes it possible to realize and appreciate that 
a very large region of the characterizing parameter S2 < s < oo is of relevance 
only to a very small physical region < g < 1 for the case where the initial 
position of the planet is inside the Schwarzschild horizon. 

For Region II, as s increases its value above those on its lower boundary 
s = si on which fc^ = 1, the value of k"^ calculated from cq.(35) decreases from 
1. The curves of constant for fc^ = 0.9,0.8,... can be easily obtained from 
eq. (35) where A and B are expressed in terms of s and e (again using MAPLE 
FSOLVE) and they arc presented in Fig. 2. However, the value of fc^ has a 
minimum value that is not in Region II. First, we show in Appendix C that 
the fc^ = 1/2 curve is given by 



The significance of this fc^ = 1/2 curve is that on it, as e ^ 1, s ^ oo, just 
like the curve for the upper boundary of Region II represented by eq.(57). For 
k^ > 0.5, the constant fc^ curves intersect the e = 1 line at some finite value of s, 
whereas for k^ < 0.5, the constant k^ curves intersect the upper boundary curve 
given by eq.(57) at points for which the values of e are less than 1. We then find 
that the minimum value of k^ in Region II is equal to 1/2-1/(2^5) = 0.276393 
which is obtained by setting e = 0, s = 1 in eqs.(32) and (35), and this value of 
k^ appears at one coordinate point only at e = and s = 1. There is no orbit 
whose k'^ is less than 0.276393 in Region II (see Fig. 2), and k"^ is thus restricted 
to the range 0.276393 < fc^ < 1. 

In Table 9, we present the coordinates (e, s) of these curves of constant fc^ 
between 0.276393 and 1. In Table 10, we present the values of q2 given by 
eq. (37) , the initial distance of the planet from the blackhole. Note that unlike 
qi for the terminating orbits in Region I whose range is finite and small, q2 can 
be infinite (for e = 1 and fc^ > 0.5). Like the terminating orbits of Region I, the 
terminating orbits of Region II can be characterized by q2 and the angle (j>2 at 
which the planet enters the center of the blackhole. If we define the "precession 
angle" Aip for the terminating orbits as in eq.(25), with k and 7 defined by 
eqs.(35) and (34), then (j)^ = K{k)/-f = A<p/2 + tt, or 



In Table 11, we present the values of (j)2- Tables 10 and 11 are to be used in 
conjunction with Table 9 that give the coordinates of the constant /c^ curves. Ex- 
amples of these terminating orbits obtained from eq.(36) are shown in Fig. 7(b)- 
(d). Again the dotted line shows the continuation of the orbit beyond 4>2- A 
planet coming from very far away, i.e. an unbounded orbit with e = 1, with 
an initial trajectory perpendicular to the line joining it to the blackhole, can 
terminate at the blackhole; the condition for this to happen is s > 0.25. Figures 
8 (a)-(c) show three unbounded orbits (e = 1) as s increases from just below 




(58) 
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to just above the critical field parameter s = 0.25. Figure 8(a) also shows an 
example of a precession angle in which the planet makes more than three revolu- 
tions around a blackhole before assuming a distance equal to its initial distance 
(which is infinity) from the blackhole. As noted after eq.(25), the actual pre- 
cession angle in this case should be more appropriately given by 2K{k)/j — Gn 
which can be obtained from the presented value of A^/tt = 4.6378 (where A<^ 
is defined by eq.(25)) and gives 0.63787r. Thus for Fig. 8(a). 0.63787r gives the 
angle between the initial incoming trajectory from very fax away at (j) = and 
the final outgoing trajectory going to infinity, i.e. 0.63787r is the polar angle 
of the direction of the outgoing trajectory going to infinity with respect to the 
a;-axis (but we have not extended the outgoing trajectory far enough to show 
the accuracy of this angle). 

We now present some useful simple expressions for the following special cases. 

(i) Special Case on the Upper Boundary of Region II Given By Eq.(57) 
We show in Appendix D that on the upper boundary of Region II given by 

eq.(57), the values of and 7 given by eqs.(35) and (34) become 

(60) 
(61) 

where the s values for 32 are given by eq.(57). 

(ii) Special Case on the Right Boundary e = 1 of Region II 

Just as for Region I, there are simple and interesting relations among fc^, s 
and 7 on the right boundary e = 1 of Region II, and they are shown in Appendix 
D. In particular, we have, on e = 1 in Region II, that for s > 1/4, 

«=^=i + i (62) 
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1 




2 


4^1/3 - 
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1/4 


7 = 


4 







or that, for 1 > A;^ > 1/2, 



and that 



^=8(fc2-l/2)' 



7=^|. (64) 

2 



Appendix D also presents a special case given by s = 1/12 that is notable. 



5 Region II' 

While we may call the entire sector s > S\ given by eq.(40) above Region I 
in Fig.2 just one region that allows only terminating orbits given by eq.(36), 
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it is useful to divide it into Regions II and IF using the curve s = S2 given by 
eq.(57). Region IF is the region of parameter space in (e, s) for which s > S2 and 
< e < 1. The heavy solid curve labeled S2 in Fig. 2 delineates the boundary 
of Region IF which separates it from Region II. Despite the apparent large 
size of Region IF, the terminating orbits here have little variety in the sense 
that the range of initial distances q2' that are given by 1/^2' = 1/3 + 4a (see 
eq.(37)); is limited (0 < g2' < 1) and the range of the angle cf)2' = K{k)/j 
at which the planet enters the blackhole is also limited. It can be shown that 
the range of <^2' is < (1)2' < O.TSQtt. An example of a terminating orbit 
obtained from eq.(36) in Region IF is shown as the solid line in Fig. 7(e); the 
dotted line shows the continuation of the orbit beyond 4>2' ■ It may be of some 
mathematical interest to note that as s — > oo in Region IF, the modulus of 
the Jacobian elliptic functions used to describe the orbits does not go to zero; 
instead fc^ ^ (2 - V3)/4 = 0.0669873, and thus P in Region IF is restricted to 
the range 0.0669873 < fc^ < 0.5. 

Tables 1-8, the orbit equations (20), (29) and (36), and the description of 
the orbits and the three regions where these orbit equations apply, complete our 
characterization of all possible planetary orbits in the Schwarzschild geometry. 

We now briefly discuss how all this may be used for the Kerr geometry 
when the spinning blackhole has a spin angular momentum per unit mass of the 
blackhole that is relatively small compared to the orbital angular momentum 
per unit mass of the planet. 



6 Kerr Geometry 

The spinning blackhole is assumed to have a spin angular momentum J given 

by [1] 



J = Mac, 



(65) 



where ac can be identified as the spin angular momentum per Tinit mass 
of the blackhole and is the quantity to be compared with h, the orbital angu- 
lar momentum per unit mass of the planet. The Kerr geometry becomes the 
Schwarzschild geometry in the limit ac/h 0. 

The worldline of a particle moving in the equatorial plane 9 — tt/2 satisfies 
the equations [1] 



1 



1 

D 



aacK 



aa 



aah 



K — 



cr 



(66) 



(67) 



where D 



ar + a^. For the equatorial trajectories of the planet in the 



Kerr geometry, the combined energy equation is 
/i^ - a^c^(fi;^ - 1) a{h-acKY 



.2 

r + 



(?a 



C^K^ - 1). 



(68) 
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Provided that a} jo? < 1 and acnjh << 1, to the first order in acn/h, it is not 
difficult to see, by comparing eqs.(66)-(68) with eqs.(3)-(6), that we can re-scale 
a to a' = a(l — 2acK,/h), s to s' = s(l — acn/h), and (pto (p' = cf)[L — b{acK/h)], 
where b is some approximation constant, such that the results we have presented 
for the orbits in the Schwarzschild geometry axe approximately applicable for 
the orbits in the Kerr geometry in terms of the scaled parameters. That is, the 
orbits in the equatorial plane and their characterization for the Schwarzschild 
and Kerr geometries are qualitatively very similar to the first order in acn/h 
except that the basic parameters s, a and 4> have to be slightly rescaled. Levin 
and Perez-Giz [5] obtained their orbits in the Kerr geometry from numerically 
integrating eqs.(66)-(68) and it would be interesting to study and examine when 
and how the planet's orbits in the Kerr geometry that they obtained can be 
related with our results with the scaled parameters, and when and how they 
begin to differ significantly from those in the Schwarzschild geometry that we 
presented in this paper. 



7 Trajectory of Light 

We now consider the deflection of light by a gravitational field. We cannot use 
the proper time t as a parameter. So we use some affine parameter a along the 
geodesic [1]. Considering motion in the equatorial plane, the geodesic equations 
give eqs.(3) and (5), and we replace the r-equation (4) by the first integral of 
the null geodesic equation, and we have [1] 



t = K, (69) 

r ■' 



= h, (71) 

where the derivative • represents d/da. Substituting eqs.(69) and (71) into 
(70) gives the 'combined' energy equation 

.%5(l-2)=cV. (72) 



„2 



Substituting dr/da = {dr/d(j)){d(j)lda) = {h/r ){dr/d(j)) and u = 1/r into 
the combined energy equation gives the differential equation for the trajectories 
of light in the presence of a gravitational field 



d'li \ ' 



The constants k and h have a physical significance through their ratio k/Zi 
as follows: Let R denote the distance of the light beam to the center of a star or 
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blackholc when the trajectory of the hght beam is such that du/dcf) = 0. R can 
either be associated with the distance of closest approach of the hght beam to 
the blackhole or with the initial distance to the blackhole of the light beam. The 
latter case is associated with light trajectories that terminate at the blackholc. 
With R so defined and letting ui = l/R, we can set c^k^ /K^ to be equal to 
u\ - au\ [8] . 

It is again convenient to consider the problem in terms of the dimensionless 
inverse distance U defined by 

a 1 

U = - = au = -. (74) 
r q 

U defined here is slightly different from the U defined by eqs.(8) and (19) 
previously. In terms of U of eq.(74), eq.(73) becomes 



Since dU /d<f) = at r = i?, one root U which we call 

Ui = ^ = awi (76) 

of the cubic equation — U"^ + c^k^o? jh? = is known, and the term 
(?K^o? jl? on the right hand side of eq.(75) can be replaced by ~U\ + 
and the other two roots of the cubic equation C/^ — [/^ — [/f + [/f = can be 
found from solving a quadratic equation. We denote the three roots of the cubic 
equation by 61,62,63. Thus writing eq.(75) as 

= C/3 - [/2 _ [/3 + [/2 (77) 

d4> ) 

the trajectory of light represented by an equation for U as a function of the 
polar angle ^ obtained from integrating eq. (77) can be characterized by a single 
parameter Ux which essentially specifies either the distance of the closest ap- 
proach or the initial distance of the light beam to the blackhole (These distances 
are scaled by the Schwarzschild radius of the blackhole) . As in our discussion of 
the planets, our references to the initial position of the light beam assume that 
the trajectory of the light beam at that initial position is perpendicular to the 
line joining that position to the star or blackhole. The range of U\ is clearly 
between and c», where J7i = means that the light beam is infinitely far away 
from the blackhole, C/i = 1 means that the light beam is at the Schwarzschild 
radius at its closest approach or its initial position, and U\ = oo means that 
the light beam is at the center of the blackhole. As we show in the following, 
the region < C/i < oc' can be appropriately divided into three sectors which 
we again call Regions I, II and IF. The similarity between the characterization 
of these three regions with that for the planetary orbits discussed in the pre- 
vious sections will bcc;onie apparent. Not surprisingly perhaps, only a single 
parameter which we choose to be I7i, is needed for the characterization of the 



23 



trajectories of a light beam in contrast to the two parameters (which we choose 
to be e and s), which we needed for the characterization of the orbits of a planet. 
The relationship between Ui and R, from eqs.(76) and (2), is 



R 



( 



GM 



Ui V c2 



Region I: < C/i < 2/3, or oo > i? > ZGMjc^ 

Here R denotes the distance of closest approach of a light beam that comes 
from a great distance. We let 



^[l-C/i + (l + 2C/i-3C/f)i/2], 



■ 2 

63 = ^[l-C/i-(l + 2C/i-3t/2)i/2], (78) 

with ei > 62 > 63, and we consider the region ei > 62 > f/ > 63, and write 
eq.(77) as 

^) ' = (61 -C/)(62 -[/)([/ -63). (79) 

Equation (79) can be integrated [7] with expressed in terms of an inverse 
sn function. After a little algebra and re-arrangement, we find the trajectory's 
equation in terms of the Jacobian elliptic functions of modulus fc to be 

1 ^ (61 - 63)62 - (62 - e3)eisn^(70, fc) 
9 (61 - 63) - (62 - e-i)8n^{^^, k) 

where 



7 



(61-63)1/^ 



62 



2 

63 



61 - 63 



(81) 



The angle of deflection Acj) can be obtained as follows. If we set g = 00 and 
also set (j) = 7r/2 + A0/2 as the incoming angle in eq.(80) (see Fig.9a for the 

special case of A(f>/2 = 45°) where A(f> denotes the total angle of deflection of 
light by the mass M, we get the following equation for determining Atj) exactly: 



sn 



7(: 



),fc 



(ei - 63)62 

(62 - 63)61 ' 



where 61, 62, 63, 7, k, are given by eqs.(78) and (81). It can also be expressed 



as 
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£i -I 

— ttH — sn~ {ip,k) 



(82) 



where 



(ei - 63)62 

(62 - 63)61 



-,1/2 



Equations (80)-(82) were first given by one of us in ref.6. Examples of these 
trajectories obtained from eq.(80) are presented in polar coordinates (g, 0) in 
Fig.9, where the blackhole is located at the origin. By setting the angles of de- 
flection A(p presented in Fig.9 to be 7r/2, tt, 37r/2, 27t, the corresponding values 
of Ui can be determined from eqs.(82) and (78) using the MAPLE FSOLVE 
program, and they are found to correspond to the distances of closest approach 
R = 4.6596GM/c2, 3.5206GM/c^, 3.2085GM/c^, 3.0902GM/c2 respectively. 
The case of R = 3.5206GM/c^ is interesting as it corresponds to the light ray 
being turned around by 180°. That the upper boundary of Region I character- 
ized by J7i = 2/3 or i? SGM/c^ is a very special case can be seen mathemati- 
cally because it results in ei = 62 = 2/3, 63 = —1/3, and hence /c^ = 1, 7 = 1/2 
and U = 2/3 = const, from eqs.(81) and (80). Physically, it results in the light 
circling the blackhole with a radius R = ZGM/c^ even though the trajectory 
has been shown to be an unstable one [1] . This known result can also be simply 
obtained from the equation of motion cPU/d(j)'^ = (3/2)C/^ — U for U = const. 
and thus U = 2/3. If one compares the size of the unstable circular photon 
orbit with the allowed limiting radii of the planetary asymptotic periodic orbits 
(2 < gmin < 2.25 or AGM/c^ < rmin < 4.5GM/c^), one can see that the radius 
of the asymptotic circular path of a planet around a blackhole is still a little 
larger than that for a photon, but not by much. 

The lower boundary of Region I characterized by ?7i = or i? = 00 gives 
ei = 1.(32 = = 0, k'^ = and 7 = 1/2, and thus gives U = or r = 00 
which is a limiting case as the light ray that is infinitely far away at its closest 
approach to the blackhole is completely undeflected. 

As in the case of the Region I particle orbits discussed in Section 3, the 
squared modulus of the elliptic functions that describe the trajectories of 
light here also covers the entire range < fc^ < 1; it varies from at the lower 
boundary to 1 at the upper boundary. 

For small Ui, the trajectory of light given by eq.(80) has been shown [6] to 
reduce to 




(83) 



and the total deflection of light to reduce to the well known result 



A(f) ~ 



AGM 
c^R ■ 



(84) 
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It can be shown from our exact result given by eq.(82) that this approximate 
expression (84) still gives an accuracy of two significant figures for Ui =0.1 or 
R = 20GM/c2. 

As Ui approaches 2/3, or as R approaches 'iGM /c? , we may let Ui = 2/3 — 5, 
where 6 = (2/3)(l — SGM/c^R) is a small positive number. From eqs.(78) 
and (81), we can express the quantities 2/7, tp and k appearing in eq.(82) 
in power series in S and find, to the first order in S, 2/7 ~ 4(1 — S/2 + ...), 
1 — 5/2+..., and k ^ 1 — 5+ .... Substituting these into eq.(82) immediately 
gives an expression for A(f) which is correct to the first order in S. If an attempt 
is made to find an expansion of sn~^(ip, k) near fc = 1, since sn~^(^j,l) = 
tanh-i V = ln[(l + ^)/(l - V')]^/^ the expansion would involve terms in In 6 
(which is a large number for small 5) and ordering the expansion terms in the 
right way can be tricky. Different forms of such expansions have been given 
and studied by various authors [9]. As we showed above and in Fig. 9, our exact 
expressions given by eqs.(80) and (82) can be used simply and directly for all 
cases in Region I. 

As Ui increases beyond 2/3 or as the distance of closest approach R of the 
light beam to the blackhole becomes smaller than 3GM/c^, the light is not just 
deflected but is absorbed by and terminates at the blackhole. It is useful to 
divide the region 2/3<[/i<ooor SGM/c^ > R> into two regions that we 
call Region II and IF that are separated by the Schwarzschild horizon, as we 
discuss below. Region II is for R from iGM /c? up to the Schwarzschild horizon, 
and Region IF is for R from the Schwarzschild horizon up to the center of the 
blackhole. 

Region II: 2/3 < JJi < 1, or 'iGMjc^ > R> 2GM/c^ 

Here R denotes the initial distance to the blackhole of the light beam which 
has initial trajectory (as (p increases from 0) perpendicular to the line joining it 
to the blackhole. 

As Ui increases beyond 2/3, Ui becomes greater than [1 — C/i + (1 + 2Ui — 
3?7i)^/^]/2, and the order of the three roots must be changed to maintain the 
inequality ei > 62 > 63. We write 



ei = Ui, 

62 = ^[l-f/i + (l + 2f/i-3i72)i/2] 



= l[l-f/i-(l + 2f/i-3[/f)i/2]. (85) 



2 



We consider the region C/ > ei > 62 > 63, and write eq.(77) as 



2 



j ={U-ei)iU-e^)iU-e3). (86) 

Equation (86) can be integrated [7] with expressed in terms of an inverse 
sn function. After some rearrangement, we find 
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1 _ ei - e2sn'^{'y(j),k) 
q cn^(7^, fc) 



(87) 



where 7 and are calculated using the same expressions given by eq.(81) 

but with ei, 62, 63 given by eq.(85). 

The expressions for 61,62,63 given by eqs.(78) and (85) coincide at fc^ = 1 
for which ei = 62 = 2/8,63 = —1/3, and both equations (80) and (87) give 
[/ = 2/3 or r = SGM/e^ independent of (j). 

Eq.(87) gives a trajectory of light which terminates at the blackhole when 
4) — 4)2 = K{k)/j. As in our discussion of the terminating orbits for the planet, 
the terminating light ray trajectories can be characterized by the angle ^2 with 
which the light beam enters the center of the blackhole. 

As C/i increases from 2/3 to 1, covers the entire range 1 > fc^ > 0; it 
decreases from 1 to 0. When Ui = 1, i.e. when the light beam grazes the 
Schwarzschild horizon, ei = 1, 62 = 63 = 0, /c^ = 0, 7 = 1/2, and we have the 
trajectory of light given by 



which gives, for (p = O.U = I or r = a, and for (j) = tt,U = ooorr = 

0, i.e. the light is absorbed at the center of the blackhole. Examples of the 
trajectories of light obtained from eqs.(87) and (88) for Ui = 5/6 = 0.83333 
{R = 2AGM/c^) and 1 (i? = 2GM/c^) in Region II arc shown as the soHd lines 
in Fig. 10 (a) and (b). The path that emerges from the center of the blackhole 
when (j) is continued beyond (^2 (shown as a dotted line in Fig. 10) again may be 
interesting if the concept of whitchole is of any physical relevance. 

When the distance R to the blackhole at (/> = is inside the Schwarzschild 
horizon, the terminating path takes on a somewhat different form as we show 
below. 

Region IV: 1< [/i < 00, or 2GM/c^ > R>0 

Here R has the same meaning as that in Region II. As J7i increases beyond 

1, i.e. when R is less than the Schwarzschild radius, 61 in eq.(85) remains real 
while 62 and 63 become complex. We now write the three roots of the cubic 
equation - U"^ - + Uf = as a, b and b given by 



1 



1 



(88) 



cos 



a 



b 



b 




Ui + i{3U^ -2Ui-iy/^] 



[/i-z(3f/f_2C/i- 1)1/2]. 



(89) 



We consider the region U > a, and write eq.(77) as 




(90) 



27 



This equation can be integrated [7] with (j) expressed in terms of an inverse 
an function. After a little algebra, we find 



1 7 + a — (7 — a)cn{^(j), k) 



(91) 



where 



and 



q 1 + cn{'-f(j), k) 

7=[C/i(3C/i-2)]i/^ (92) 



1 3j7i-l 1 3a-l 

2 4VC/i(3[/i -2) 2 472 ■ 

Equation (91) gives the trajectory of light when R is inside the Schwarzschild 
horizon and it terminates at the blax;khole when (j) = 4>2' = '2K{k)/j, where k 
and 7 are given by eqs.(93) and (92). On the boundary with Region II where 
Ui = 1, and fc2 = 0, 7 = 1 from eqs.(93) and (92), eq.(91) becomes eq.(88) and 
thus there is no discontinuity in the orbit as it makes a transition from Region 
II to Region IF across Ui = 1. 

We note that as in the case of Region IF for the planetary orbits. Region 
IF for light trajectories covers a semi-infinite range of the parameter character- 
izing it (1 < t/i < 00) but is of relevance only to a very small physical region 
2GM/ > > for the initial position of a light beam inside the Schwarzschild 
horizon. The terminating orbits of light rays in Region IF are also of very little 
variety as 4)21 is restricted to a limited range of < 4)21 < tt. An example 
of a terminating trajectory obtained from eq.(91) is shown as the solid line in 
Fig.lO(c) for U\ = 10 (i? = 0.2GM/c^); the dotted line again represents a tra- 
jectory of light coming out from the center of the blackhole as is continued 
beyond ^2'- ^ may be of some mathematical interest to note that as t/i ^ cxd or 
R^Q, the squared modulus of the Jacobian elliptic functions used to describe 
the trajectories k"^ approaches a value (2 — \/3)/4 = 0.0669873 that is the same 
as that given in Section 5 for the case of Region IF for the planetary orbits. Thus 
the squared modulus of the elliptic functions that describe the terminating light 
trajectories in Region IF is restricted to a very small range < /c^ < 0.0669873 
even as Region IF consists of a very large interval 1 < C/i < 00. 



8 Summary 

We have presented exact analytic expressions given by eqs.(20), (29) and (36) 
for the planetary orbits in the Schwarzschild geometry. The equations relate the 
distance r of the planet from the star or blackhole to the polar angle (f) and are 
described explicitly by Jacobian elliptic functions of modulus k. Equation (20) 
gives periodic orbits that describe a planet precessing around a star or blackhole, 
while eqs.(29) and (36) give terminating orbits that describe a planet phmging 
into the center of a blackhole. One of the most important aspects of our analysis 
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is the construction of a map with coordinates (e, s) that we use to view all 
possible orbits in their entirety, where the two dimensionless parameters e and 
s are defined by eqs.(ll) and (12) which we call the energy and field parameters 
respectively. For < e < 1, we show that there arc three regions which we 
call Regions I (0 < s < Si), II (si < s < S2) and IF (s2 < s < 00) where 
these orbits are applicable, and where si and S2 that depend on e are given by 
eqs.(40) and (57) respectively (Fig. 2). Region I has periodic and terminating 
orbits given by eqs.(20) and (29). Regions II and IF have terminating orbits 
only given by eq.(36). We have divided Region I into grids that consist of lines 
of constant precession angle < A0 < 00 given by cq.(25) and lines of constant 
true eccentricity < e < 1 defined by eq.(28) (Fig. 3); the lines of constant At/) 
are obtained from solving eqs.(38) and (99), and those of constant e from solving 
eqs.(38) and (100). These grids make the identification of all possible periodic 
orbits, including the unbounded and asymptotic ones, convenient and precise. 
Numerous numerical results for orbits in Region I are presented in Tables 1-6, 
and examples of preccssing orbits, including the unbounded ones, are shown in 
Figs. 4 and 5. Among the interesting results, for example. Table 6 for At/) = 27r, 
e = 1, s = 0.248804 and Fig. 4 f3 show that a planet coming from infinity at 
zero polar angle, makes a complete loop about the blackholc, and returns to 
infinity with a polar angle that approaches 27r. Figure 3, or a more refined 
version of it that can be constructed using the expressions for and e that 
we presented, can be fruitfully used with the experimental observation data. 
As for the terminating orbits, our orbit equations assume the planet's initial 
trajectory to be perpendicular to the line joining it to the star or blackhole. The 
terminating orbits of Region I require the planet to be at a distance iji < 2 from 
the star or blackhole, and those of Region IF requires ^2' < 1) i-e- they require 
the planet to be initially very close to the blackhole or within the Schwarzschild 
horizon. The terminating orbits of Region II. on the other hand, arc more 
interesting as the planet can be initially at a distance 1 < 92 < 00. We have 
shown how the periodic orbits of Region I become the asymptotic periodic orbits 
as s si, and then become terminating orbits as s becomes greater than si. 
It is interesting to note that a planet initially at infinity with e = 1 and in a 
trajectory perpendicular to the line joining it to a blackhole can be absorbed 
by the blackhole if s is greater than 0.25. 

We have also presented exact analytic expressions given by eqs.(80), (87) 
and (91) for the trajectories of light in the presence of a star or blackhole 
depending on the value of one parameter Ui that has a range which can be 
divided into three regions: Regions I (0 < [/i < 2/3), II (2/3 < Ui < 1) and 
IF (1 < J/i < 00), where Ui is defined by eq.(76). In Region I, the deflection 
of light can range from small angles to going continuously around the star or 
blackhole in a circle. In Regions II and IF, light is absorbed into the center 
of the blackhole. Among the interesting results, a deflection of light by 180° 
requires a distance of closest approach R to the blackhole equal to 3.5206GM/c^ 
(Fig. 9b), and for R < SGM/c^ , light will be absorbed by the blackhole. 

We have thus presented a complete map that can help to identify character- 
istics of stars and blackholes (that are not spinning too fast) from the observed 



29 



characteristics of objects or light beams that are affected by them. 
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10 Appendix A Relation Among s, e and /c^ 

In this Appendix, we derive the relation among s,e and given by eqs.(38) 
and (39). Substituting eq.(23) into eq.(22) and after a little algebra, we find 

tan3 = ^, (94) 

and hence we find 



sin - = 



cos 



V3P 

2Vl - fc2 + ' 
2-fc^ 

2Vl - fc2 + fc^ ' 



(95) 
(96) 



We then find 



'e 47r 
3 + T 



COS 



d 27r 



-1 + 2P 



2Vl-fc2 + fc4' 

-1 - A-^ 
2Vl - A;2 + A;4 ■ 



(97) 
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Equations (96)-(98) and (23) allow 61,62,63 to be expressed in terms of 52 

and A:^ as 



61 



62 = 



2-P 
12 VI - P + fc4 ' 
fg^ -1 + 2k^ 
12 VI- fc2 + A;4' 



Substituting eq.(96) into the relation cos^ = — 3cos(0/3) +4cos^(0/3), and 
substituting the result into eq.(24) gives the relation among s,e and k"^ given 
by eqs.(38) and (39). For fc^ = 1, wc get cos 9 = -1 from eqs.(38) and (39), and 
from eq.(98) we get ei = 62 = ^Jg^J^ = —(1/2)63. 

From eqs.(25), (21), (98) and (10), we find an expression for the precession 
angle given by eq.(25) in terms of k and s to be 

A.^ = ( i_i2,2 ) - 2^^. (99) 

11 Appendix B The Energy Parameter e and 
the True Eccentricity e in Region I 

In this Appendix, wc show the relation between the energy parameter e and 
the true eccentricity e. The energy parameter e in all three regions I, II and IF 
is defined by eq.(ll). The general or true eccentricity e is defined by eq.(28). 
From cq.(28), (26) and (27) and using the expressions (98) of ei, 62, 63 in terms 
of k^ given in Appendix A, we find that £ can be expressed as 

6(62 - 63) 



1 + 6(e2 + 63) 

or 



3^'^^^ (100) 



2Vl - fc2 + fc4 _ (2 _ fc2yi_l2s2' 

For small and fc^, we find that 



1 - (1 - 6s2) 4s2 

because fc^ ~ 4532 j^j. ^^^^ small A;^ and [6]. We thus confirm the iden- 
tification of e defined in eq.(ll) with the eccentricity of the orbit in Newtonian 
mechanics. As we pointed out in the text, £ is generally not equal to e. How- 
ever, as we show below, e = e exactly when e = 1 and in this case the orbits 
are unbounded. 
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For the possibility of unbounded orbits in Region I, we set the initial q = oo 
at = in eq.(20), and get 



or 



^+463 = 0, (102) 



^3 = -^. (103) 



Using eqs.(23), (10) and (97) that give 63 in terms of s and fc^, and after a 
Uttle algebra, we find the simple equation that relates s to fc^ on e = 1 to be 
given by eq.(51). 

This simple equation (51) between s and for e = 1 can be used for any 
< A;^ < 1. For example, we find that for fc^ = 1, s = 1/4 = 0.25, and for 
A;2 = 1/2, s = 1/(3^2) = 0.235702. For e = 1, and for small fc^ we have 

~ A;^/4 which is a special case of k'^ ~ 4es^ that is valid more generally for 
< e < 1. 

In addition, we find 32 = 1/12 - = 1/12 - k^/[A{l + k'^f] and therefore 
and from the expression for ei and 62 given by eq.(98), we find 



2-fc^ 
~ 12(1 + fc2)' 
-l + 2fc2 

= mrwy ^'''^ 

Tims from the expressions (27), (30), and (21) for gmin, Qi, and 7, we find 
that when e = 1, they have the simple expressions given by eqs.(53), (54) and 
(52). 

Also, substituting eq.(48) into eq.(lOO) shows that e = 1 when e = 1, i.e. e 
and e coincide at e = 1. 



12 Appendix C Some Simple Relations for the 
Special Case of /c^ = 1/2 

It is known that in elliptic functions, the squared modulus /c^ = 1/2 is a special 
value for which many simple relations arise. We first consider the case of = 
1/2 in Region I. We note that substituting fc^ = 1/2 into eq.(22) gives a relation 
ei + 63 = 2e2, and substituting the expressions of 61,62,63 from eq.(23) into 
this relation gives 9/3 = n/6, or 6 — n/2. Thus cos 9 — which results in 

93 = (106) 



32 



from eq.(24), which in turn gives a simple relationship 



For example, we have s = ^J{3-V3)/18 = 0.265408 for e = 0, and 

s = l/(3v^) = 0.235702 for e = 1 (using L'Hospital rule). These two values 
represent the two terminal coordinates of the constant fc^ = 1/2 line in Region 
I (see Fig.l). The other solution of eq.(106) is eq.(58) which is applicable for 
Region II as we shall show later in this Appendix. 
We also find from eq.(23) that 



ei = -63 = 



and 



/92 

2 ' 

62 = 0, (108) 



(109) 

and the orbit equation (20) becomes 
Q <J 

where the s value for 92 in this case is given by cq.(107) for Region I. The 
precession angle can be found from eq.(25) and from K{1/^) = 1.85407. 
It is given by 

Acj) _ 1.18034 _ 2 



From eqs.(26)-(28) and (108), we find 

6(62-63) ^ 3Vfl^ 
1 + 6(62 + 63) 1-3^' ^ ' 

which can be inverted and solved for s in terms of e, giving 

2 3 + 6e-£2 
^=36(TW 
Substituting eq.(lll) into eq.(106) gives e in terms of s as 

12(l + £)^(l-£)(l + 3g) 
^ (3 + 6£-£2)2 • ^^^^> 

Since for 6 = 0, s = ^J{3-VS)/18 when fc^ = 1/2 as we showed above, sub- 
stituting this s value into eq.(106) gives e = (2/\/-3 + 2^3 - 1)"^ = 0.516588. 
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For e = 1, s = l/(3\/2) when fc^ = 1/2, and substituting this s value into 
eq.(106) gives £ = 1 as it should. Thus eqs.(lll) and (112) are the parametric 
equations for the line of constant fc^ = 1/2 which can be used instead of eq.(107) 
as £ takes the values between 0.516588 and 1. 

Also, substituting 62 = from eq.(108) into eq.(27) gives 

gmin = 3 (113) 

independent of e for fc^ = 1/2, as shown in Table 3. 

We now consider = 1/2 in Region II. From eqs.(35) and (33), we have 

a = A+B = 0. (114) 

From cq.(32), and from A = —B, and = —B^, wc arrive again at cq.(106) 
with the same (73 given by eq.(lO). This explains why we stated after eq.(107) 
that the other solution of eq.(106) given by eq.(58) gives the relation between s 
and e for Region II. The two equations (108) and (58) that give simple relations 
between s and e in two different regions and that arise as two different solutions 
of the same equation (106) show a rather remarkable symmetry exhibited by 
the special case = 1/2. 

It also follows from eqs.(32) and (34) that 



lilts' -1) 



1/4 



(116) 



and that the orbit equation (36) becomes 



q 3 l + cn(27(/.,l/\/2)' ^ ' 

where the s value for the above equations is given by eq.(58) for < e < 1. 
Since a = 0, the initial distance q2 of the planet from the blackhole is 92 = 3 
from eq.(37), independent of e, as shown in Table 10. 

13 Appendix D The Boundaries of Region II 

On the upper boundary 53 = 1/(1 ~ e^) of Region II, the planet starts from 
the Schwarzschild horizon given by g = 1, which implies that for = 0, 1/q = 
1/3 + 4a from eq.(36), or 



I + 4a. (118) 



Hence 
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a = A+B=^. (119) 



Since from eqs.(32) and (17), 



wc can conclude from eqs.(35) and (34) that k'^ and 7 are given on the 
boundary = 1/(1 ~ e^) of Region II by eqs.(60) and (61). 

In particular, for e = 0, s = 1, we find 7 = (5/16)1/'' = 0.747674, P = 
(1 - l/\/5)/2 = 0.276393, and this is the minimum value of k'^ in Region II. 

We now consider the right boundary e = 1 of Region II. Consider the un- 
bounded orbit of a planet coming from infinity that requires, from eq.(36), that 



3 

or 



^+4a = 0, (121) 



a=-^. (122) 



Prom 



A + B = a = -^, (123) 



12 
and 

AB=^, (124) 

and from cqs.(34) and (35), wc find cqs.(62) and (64) that give and 7 in 
terms of s on the boundary e = 1 of Region II. 

Finally, there is a special case of 172 = or = 1/12 (s = l/2\/3 = 
0.288675135) in Region II that is somewhat interesting. If (72 of eq.(lO) is 
equal to zero, then A = 27s'3 from eq.(17), and from eqs.(32)-(35), we find 
A = (2^3)1/3/2, B = 0, a = A, J = (3^^)1/4^ and noting that A < for 

= 1/12 and < e < 1, we have 

fc2 ^ 1 _ V3A ^ 1 %/3 ^ 0.933012702 

2 4 I A I 2 4 

which is independent of e, i.e. the constant fc^ = (2 + -s/S) /4 curve in Region 
II just above the boundary curve of Regions I and II in Fig.2 is a horizontal 
line. Thus the terminating orbits represented by eq.(36) for any point along this 
horizontal line are represented by elliptic functions of the same squared modulus 
given above. It should be remembered, however, that the orbits and the initial 
positions of the planet (the initial trajectory of which is perpendicular to the 



35 



line joining it to the blackhole) depend also on 7 and a that are dependent on 
the value of e. Thus the orbit for e = 1, for example, that gives a = —1/12 
(and 7 = 1/(2 • 3^/'') = 0.379917843) is a terminating orbit for a planet that is 
initially infinitely far away from a blackhole. 
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Table 1 : Values of s for various values of e and in Region I 





s 


e = 0.0 


e = 0. 1 


(? = 0.2 


e = 0.3 


e - 0.4 


e = 0.5 


e = 0.6 


e = 0.7 


e = 0.8 


c = 0.9 


e = 1.0 


e 


= 0.001 


0.0554787 


0.0454170 


0.0346617 


0.0286392 


0.0248892 


0.0222952 


0.0203686 


0.0188667 


0.0176539 


0.0166480 


0.0157963 


e 


-a 01 


0.114809 


0.109191 


0.0969713 


0.0851063 


0.0757980 


0.0686999 


0.0631762 


0.0587575 


0.0551328 


0.0520954 


0.0495050 


e 


= 01 


0.210213 


0.208385 


0.203281 


0.195887 


0.187374 


0.178683 


0.170376 


0.162703 


0.155729 


0.149428 


0.143740 


e 


= 02 


0.238703 


0.237612 


0.234488 


0.229734 


0.223875 


0.217424 


0.210787 


0.204236 


0.197932 


0.191956 


0.186339 


e 


-03 


0.252575 


0.251809 


0.249595 


0.246160 


0.241814 


0.236873 


0.231612 


0.226239 


0.220900 


0.215689 


0.210663 


k' 


-04 


0.260533 


0.259944 


0.258236 


0.255562 


0.252131 


0.248163 


0.243857 


0.239372 


0.234829 


0.23031 1 


0.225877 


e 


= 05 


0.265408 


0.264926 


0.263523 


0.261314 


0.258458 


0.255122 


0.251460 


0.247600 


0.243642 


0.239658 


0.235702 


e 


= 06 


0.268462 


0.268045 


0.266831 


0.264913 


0.262422 


0.259494 


0.256258 


0.252821 


0.249269 


0.245666 


0.242061 


e 


= 07 


0.270350 


0.269973 


0.268875 


0.267137 


0.264873 


0.262203 


0.259239 


0.256076 


0.252791 


0.249443 


0.246076 


e 


-OS 


0.271452 


0.271099 


0.270069 


0.268436 


0.266305 


0.263787 


0.260985 


0.257986 


0.254863 


0.251671 


0.248452 




= 09 


0.272006 


0.271665 


0.270668 


0.269088 


0.267025 


0.264583 


0.261863 


0.258948 


0.255908 


0.252796 


0.249653 


e 


= 1.0 


0.272166 


0.271828 


0.270840 


0.269276 


0.267232 


0.264812 


0.262116 


0.259225 


0.256209 


0.253120 


0.250000 



Table 2; Values of for various values of e and lO- in Region I 



e-OO 


e-O.l 


e.0.2 


e-0.3 


e-OA 


6 = 0.5 


6-0.6 


e = 0.7 


6 = 0.8 


€-0.9 


e - 1.0 


174.74 


273.42 


522.17 


871.99 


1345.9 


2012.0 


3012.8 


4681.4 


8018.4 


18025 


00 


43.673 


49.688 


68.266 


99.742 


146.78 


212.38 


313.52 


482.96 


822.60 


1842.4 


00 


14.306 


14.851 


16.550 


19.608 


24.429 


31.802 


43.389 


63.171 


103.21 


223.98 


00 


1 1 .377 


1 1 .697 


12.693 


14.481 


17.300 


21 .627 


28.458 


40.168 


63.938 


135.75 


00 


10.282 


10.534 


11.314 


12.713 


14.914 


18.291 


23.625 


32.774 


51.357 


107.53 


00 


9.7280 


9.9479 


10.630 


1 1 .849 


13.766 


16.704 


21.342 


29.298 


45.460 


94.318 


00 


9.4118 


9.6145 


10.243 


11.365 


13.128 


15.828 


20.089 


27.397 


42.239 


87.109 


00 




9.4148 


10.012 


11.078 


12.751 


15.314 


19.356 


26.285 


40.359 


82.903 


00 


9.1078 


9.2946 


9.8733 


10.906 


12.527 


15.008 


18.920 


25.627 


39.246 


80.413 


00 


9.0421 


9.2256 


9.7938 


10.808 


12.399 


14.833 


18.672 


25.252 


38.612 


78.997 


00 


9.0094 


9.1912 


9.7542 


10.759 


12.335 


14.747 


18.549 


25.066 


38.299 


78.296 


00 


9.0000 


9.1814 


9.7429 


10.745 


12.317 


14.722 


18.514 


25.013 


38.209 


78.095 


00 



Table 3: Values of ^mn for various values of e and in Region I 



„i„ (- = 0.0 f = 0. 1 f = 0.2 f = 0.3 f = 0.4 e = 0.5 r = 0.6 e = 0.7 r = 0.8 e = 0.9 e = 1.0 

0.001 149.15 215.27 343.84 466.77 574.79 669.11 751.92 825.10 890.20 948.48 1001.0 

OOl 31.135 33.980 41.469 50.945 60.399 69.110 76.948 83.963 90.248 95.899 101.00 

= ai 7.0549 7.1489 7.4154 7.8125 8.2864 8.7891 9.2876 9.7635 10.209 10.621 11.000 

= 02 4.7480 4.7753 4.8530 4.9703 5.1133 5.2691 5.4278 5.5830 5.7311 5.8703 6.0000 

= 03 3.8359 3.8465 3.8765 3.9220 3.9777 4.0389 4.1018 4.1639 4.2237 4.2803 4.3333 

= 04 3.3289 3.3325 3.3428 3.3583 3.3773 3.3983 3.4199 3.4413 3.4619 3.4815 3.5000 

= 05 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 3.0000 

= 06 2.7667 2.7645 2.7585 2.7494 2.7382 2.7260 2.7134 2.7009 2.6889 2.6774 2.6667 

= 07 2.5911 2.5877 2.5778 2.5629 2.5446 2.5247 2.5042 2.4840 2.4645 2.4460 2.4286 

= 08 2.4535 2.4491 2.4366 2.4178 2.3949 2.3698 2.3442 2.3189 2.2946 2.2716 2.2500 

= 09 2.3422 2.3371 2.3228 2.3012 2.2751 2.2465 2.2174 2.1887 2.1613 2.1354 2.1111 

= 1.0 2.2500 2.2445 2.2288 2.2052 2.1767 2.1458 2.1142 2.0833 2.0538 2.0259 2.0000 
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Table 4: Values of h^ln for various values of e and in Region I 







e = 0.0 




-0.001 


0.018906 


k' 


-0.01 


0.088019 


k 


= 0.1 


0.4221 1 


k 


= 0.2 


0.69762 


k 


-0.3 


0.95650 


k 


= 0.4 


1.2200 


k 


= 0.5 


1.5029 


k 


-0.6 


1.8226 


k 


= 0.7 


2.2072 


k 


= 0.8 


2.7169 


k 


= 0.9 


3.5401 


k 


= 1.0 


00 



e-0.1 6 = 0.2 

0.012572 0.0072747 

0.078760 0.060819 

0.41041 0.37969 

0.68449 0.64888 

0.94206 0.90245 

1.2043 1.1608 

1.4858 1.4383 

1.8038 1.7519 

2.1866 2.1293 

2.6938 2.6295 

3.5129 3.4374 

OO QO 



c = 0.3 c-0.4 

0.0049522 0.0037347 

0.046033 0.036087 

0.33964 0.29902 

0.59968 0.54604 

0.84652 0.78377 

1.0987 1.0281 

1.3701 1.2917 

1.6770 1.5905 

2.0465 1.9507 

2.5365 2.4286 

3.3281 3.2011 

QO QO 



c = 0.5 e-0.6 

0.0029941 0.0024975 

0.029413 0.024739 

0.26259 0.23177 

0.49423 0.44737 

0.72123 0.66290 

0.95649 0.88857 

1.2116 1.1348 

1.5015 1.4157 

1.8517 1.7559 

2.3170 2.2087 

3.0696 2.9419 

QO QO 



e-0.7 e-0.8 

0.0021419 0.0018748 

0.021316 0.018713 

0.20625 0.18522 

0.40639 0.37108 

0.61042 0.56410 

0.82648 0.77085 

1.0638 0.99972 

1.3360 1.2635 

1.6666 1.5851 

2.1075 2.0150 

2.8224 2.7130 

QO 00 



e-0.9 e = 1.0 

0.0016669 0.0015004 

0.016671 0.015029 

0.16780 0.15323 

0.34079 0.31478 

0.52353 0.48808 

0.72150 0.67786 

0.94234 0.89123 

1.1982 1.1398 

1.5114 1.4453 

1.9312 1.8558 

2.6139 2.5247 

00 00 



Table 5: Values of e for various values of e and fe^ in Region I 





= 0.001 


0.079005 


k^ 


= 0.01 


0.16760 




= 0.1 


0.33947 




= 0.2 


0.41109 




= 0.3 


0.45659 




= 0.4 


0.49009 




= 0.5 


0.51659 




- 0.6 


0.53845 




= 0.7 


0.55703 




= 0.8 


0.57314 




= 0.9 


0.58733 


e 


= 1.0 


0.60000 



e = 0. 1 
0.11900 
0.18774 
0.35009 
0.42021 
0.46503 
0.49813 
0.52436 
0.54603 
0.56445 
0.58044 
0.59454 
0.60713 



p = 0. 2 
0.20592 
0.24420 
0.38115 
0.44683 
0.48963 
0.52152 
0.54692 
0.56798 
0.58594 
0.60155 
0.61534 
0.62766 



e = 0. 3 
0.30268 
0.32383 
0.43017 
0.48894 
0.52845 
0.55832 
0.58231 
0.60232 
0.61944 
0.63437 
0.64759 
0.65942 



f = 0.4 
0.40148 
0.41413 
0.49343 
0.54373 
0.57889 
0.60598 
0.62797 
0.64644 
0.66233 
0.67623 
0.68857 
0.69963 



e = 0. 5 
0.50088 
0.50896 
0.56695 
0.60819 
0.63825 
0.66190 
0.68133 
0.59778 
0.71200 
0.72449 
0.73560 
0.74558 



p = 0. 6 
0.60054 
0.60587 
0.64738 
0.67964 
0.70412 
0.72378 
0.74014 
0.75410 
0.76623 
0.77692 
0.78645 
0.79502 



e = 0.7 
0.70032 
0.70380 
0.73227 
0.75594 
0.77454 
0.78978 
0.80261 
0.81364 
0.82327 
0.83179 
0.83939 
0.84623 



e = 0. 8 
0.80015 
0.80227 
0.81998 
0.83548 
0.84802 
0.85847 
0.86737 
0.87508 
0.88183 
0.88781 
0.89316 
0.89798 



e = 0.9 
0.90002 
0.90105 
0.90946 
0.91710 
0.92344 
0.92880 
0.93341 
0.93743 
0.94096 
0.94410 
0.94690 
0.94943 



Table 6: Values of s for constant values of in Region I 





s 


e = 0. 


e = 0. 1 


e = 0. 2 


(• = 0.3 


e = 0. 4 




= nils 


0.0929975 


0.0929922 


0.0929763 


0.0929498 


0.0929129 




= nl6 


0.150971 


0.150946 


0.150871 


0.150747 


0.150574 




= 1/3 


0.195246 


0.195185 


0.195000 


0.194694 


0.194273 


A^ 


= nil 


0.220477 


0.220377 


0.220080 


0.219591 


0.218920 


A,^ 


= K 


0.254214 


0.254018 


0.253437 


0.252492 


0.251216 


A0 = 


inl2 


0.265371 


0.265111 


0.264346 


0.263113 


0.261468 


A^ 


= 2k 


0.269502 


0.269206 


0.268334 


0.266938 


0.265091 


A^ 


= 00 


0.272166 


0.271828 


0.270840 


0.269276 


0.267232 



p = 0. 5 


e = 0.6 


e = 0. 7 


(' = 0. 8 


e = 0. 9 


e = 1.0 


0.0928655 


0.0928078 


0.0927400 


0.0926622 


0.0925745 


0.0924773 


0.150353 


0.150087 


0.149776 


0.149424 


0.149031 


0.148601 


0.193741 


0.193104 


0.192372 


0.191551 


0.190651 


0.189680 


0.218080 


0.217085 


0.215951 


0.214696 


0.213336 


0.211888 


0.249650 


0.247838 


0.245823 


0.243650 


0.241354 


0.238971 


0.259478 


0.257210 


0.254729 


0.252091 


0.249347 


0.246537 


0.262877 


0.260379 


0.257671 


0.254819 


0.251875 


0.248804 


0.264812 


0.262116 


0.259225 


0.256209 


0.253120 


0.250000 
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Table 7: Values of s for constant values of e in Region I 



s «-0.00 e-0.01 e-O.m e - 0.06 c = 0.08 c = 0.10 c = 0.12 c-0.14 6-0.16 e-0.18 « - 0.20 

e = 0.2 0.135153 0.134536 0.132663 0.129457 0.124776 0.118388 0.109916 0.0987083 0.0835085 0.0611485 

s e = 0.00 « = 0.04 « = 0.08 c = 0.12 c = 0.16 c = 0.20 c = 0.24 e = 0.28 « = 0.32 £-0.36 e = 0.40 

e = 0.4 0.234806 0.234069 0.231803 0.227832 0.221831 0.213264 0.201255 0.184309 0.159570 0.120022 

s « = 0.00 « = 0.06 « = 0.12 c = 0.18 c = 0.24 c = 0.30 c = 0.36 e - 0.42 6-0.48 e = 0.54 « = 0.60 

e = 0.6 0.272166 0.272037 0.271572 0.270524 0.268439 0.264552 0.257565 0.245130 0.222468 0.177299 

s « = 0.6I c = 0.62 e-0.64 c = 0.66 c = 0.68 c = 0.70 c = 0.72 e - 0.74 e - 0.16 e - 0.7S « = 0.80 

£-0.8 0.261834 0.261450 0.259987 0.257309 0.252931 0.246121 0.235710 0.219699 0.194219 0.149735 



Table 8: Values of qi for various values of e and k'^ In Region I {terminating orbits) 





91 


«=0.0 


e = 0.1 


p. 


= 0.001 


1.0126 


1 .0084 




= 0.01 


1 .0582 


1 .0521 




= 0.1 


1 .2685 


1.2614 




= 0.2 


1 .4255 


1.4182 




= 0.3 


1 .5575 


1 .5502 




= 0.4 


1 .6756 


1 .6683 




= 0.5 


1 .7844 


1 .7773 


k^ 


= 0.6 


1 .8864 


1 .8795 


k^ 


= 0.7 


1 .9831 


1 .9764 




= 0.8 


2.0754 


2.0691 




= 0.9 


2.1642 


2.1583 




= 1.0 


2.2500 


2.2445 



« = 0.2 c = 0.3 c = 0.4 

1.0048 1.0033 1.0025 

1.0403 1.0306 1.0240 

1.2427 1.2180 1.1928 

1.3983 1.3703 1.3394 

1.5299 1.5007 1.4672 

1.6481 1.6186 1.5841 

1.7574 1.7281 1.6935 

1.8601 1.8315 1.7973 

1.9578 1.9301 1.8968 

2.0513 2.0248 1.9929 

2.1415 2.1164 2.0860 

2.2288 2.2052 2.1767 



c = 0.5 c = 0.6 c = 0.7 

1.0020 1.0017 1.0014 

1.0196 1.0165 1.0142 

1.1699 1.1504 1.1341 

1.3089 1.2810 1.2563 

1.4332 1.4008 1.3711 

1.5483 1.5135 1.4809 

1.6570 1.6210 1.5869 

1.7610 1.7248 1.6900 

1.8612 1.8254 1.7908 

1.9584 1.9236 1.8898 

2.0532 2.0198 1.9872 

2.1458 2.1142 2.0833 



c = 0.8 e = 0.9 e- 1.0 

1.0012 1.0011 1.0010 

1.0124 1.0111 1.0100 

1.1206 1.1094 1.1000 

1.2348 1.2161 1.2000 

1.3445 1.3209 1.3000 

1.4511 1.4242 1.4000 

1.5552 1.5263 1.5000 

1.6575 1.6275 1.6000 

1.7582 1.7279 1.7000 

1.8577 1.8278 1.8000 

1.9562 1.9271 1.9000 

2.0538 2.0259 2.0000 



Table 9: Values of s for various values of e and k^ in Region II 





s 


e = 0. 


e = 0. 1 


e = 0. 2 


e = 0. 3 


e = 0.4 


k- 


= 1.0 


0.272166 


0.271828 


0.270840 


0.269276 


0.267232 


k- 


= 0.9 


0.297739 


0.297917 


0.298442 


0.299293 


0.300443 


k' 


= 0.8 


0.329945 


0.330718 


0.333024 


0.336843 


0.342169 


k- 


= 0.7 


0.371926 


0.373406 


0.377868 


0.385401 


0.396232 


k- 


= 0.6 


0.429234 


0.431592 


0.438767 


0.451100 


0.469328 


k- 


= 0.5 


0.512730 


0.516250 


0.527046 


0.545908 


0.574478 


k- 


= 0.4 


0.646974 


0.652192 


0.668307 


0.696858 


0.741019 




= 0.3 


0.901890 


0.910054 


0.935412 


0.980848 


1.05231 



t' = 0.5 e = 0.b c = 0.7 c = 0.8 c = 0.9 e = 1.0 

0.264812 0.262116 0.259225 0.256209 0.253120 0.250000 

0.301865 0.303537 0.305444 0.307575 0.309927 0.312500 

0.349050 0.357613 0.368098 0.380913 0.396727 0.416667 

0.410826 0.430060 0.455531 0.490285 0.540844 0.625000 

0.494871 0.530424 0.581374 0.659978 0.802847 1.25000 

0.615920 0.676462 0.769603 0.930895 1.30267 

0.806949 0.907116 1.07001 1.37756 2.21530 
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Table 10: Values of qi for various values of e and in Region II 





32 


e-0.0 


e-0.1 


e-0.2 


6 = 0.3 


6 = 0.4 


6 = 0.5 


6 = 0.6 


6 = 0.7 


6 = 0.8 


6 = 0.9 


6 = 1.0 




- 1.0 


9.0000 


9.1814 


9.7429 


10.745 


12.317 


14.722 


18.514 


25.013 


38.209 


78.095 


00 


k' 


-0.9 


7.6760 


7.7965 


8.1685 


8.8291 


9.8596 


11.426 


13.833 


18.073 


26.550 


52.109 


00 


e 


= 0.8 


6.4074 


6.4806 


6.7060 


7.1041 


7.7200 


8.6471 


10.084 


12.507 


17.353 


31.823 


00 


e 


-0.7 


5.1997 


5.2381 


5.3557 


5.5620 


5.8780 


6.3467 


7.0597 


8.2346 


10.522 


17.148 


00 


e 


-0.6 


4.0606 


4.0750 


4.1191 


4.1959 


4.3120 


4.4809 


4.7310 


5.1276 


5.8592 


7.8120 


00 


e 


= 0.5 


3.0000 


3.0000 


3.0000 


3.0000 


3.0000 


3.0000 


3.0000 


3.0000 


3.0000 


3.0000 




e 


= 0.4 


2.0324 


2.0259 


2.0059 


1.9717 


1.9217 


1.8523 


1 .7574 


1 .6241 


1 .4235 


1 .0692 




e 


-0.3 


1.1805 


1.1733 


1.1516 


1.1148 


1.0618 















Table 1 1 : Values of i^iln for various values of e and in Region II 



( 


hl^ 


6 = 0.0 


6 = 0.1 


6 = 0.2 


6 = 0.3 


6 = 0.4 


6 = 0.5 


6 = 0.6 


6 = 0.7 


6 = 0.8 


6 = 0.9 


6 = 1.0 




= 1.0 


oo 


00 


oo 




CO 


oo 


oo 


00 


00 


oo 


00 




= 0.9 


2.6599 


2.6469 


2.6099 


2.5551 


2.4890 


2.4176 


2.3448 


2.2733 


2.2043 


2.1385 


2.0761 




= 0.8 


2.1586 


2.1480 


2.1175 


2.0711 


2.0132 


1.9479 


1 .8781 


1.8055 


1 .7309 


1.6541 


1.5742 




= 0.7 


1.8169 


1 .8079 


1.7819 


1.7412 


1.6889 


1.6274 


1 .5584 


1.4822 


1 .3977 


1.3010 


1.1817 




= 0.6 


1.5360 


1 .5283 


1.5060 


1.4703 


1.4231 


1.3654 


1 .2977 


1.2186 


1.1237 


1.0002 


0.78496 




= 0.5 


1.2822 


1 .2758 


1.2569 


1.2262 


1.1843 


1.1315 


1 .0671 


0.98815 


0.88720 


0.74053 






= 0.4 


1.0373 


1 .0321 


1.0166 


0.99105 


0.95540 


0.90923 


0.85112 


0.77757 


0.68006 


0.53257 






= 0.3 


0.78606 


0.78213 


0.77030 


0.75047 


0.72230 
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Fig. 1. Region I plots of = 0.001,0.01,0.1,0.3,0.5,0.7,1.0 




Fig. 2. Region II plots o\k^ = 1.0,0.9,0.8,0.7,0.6,0.5,0.4,0.3 
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Fig. 4. Region I: Periodic(al-fl,a2-f2), unbounded(a3-g3), asymptotic periodic orbits(gl-g3) for A(p equal to 
(a) 7r/6, (b) tt/S, (c) 7r/2, (d) n, (e) 3n/2, (f) 2n, (g) oo and for e equal to (1) 0, (2) 0.5, (3) 1. 
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Fig. 6. Terminating orbits in region I 

(a) = 0.01, e = 0, s = 0.114809, ^ = 1.0440, = 1.0582 

(b) = 0.1, e = 0.5, s = 0.178683, ^ = 1.1313, = 1.1699 

(c) k^ = 0.9, e = 0.9, s = 0.252796, ^ = 2.3070, q^ = 1.9271 

(d) k^ = 1, e = 1, s = 0.25, q)^ = oo, <?i = 2 
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Fig. 7. Terminating orbits in Region II (a-d) and Region IF (e) 

(a) = 1, e = 1, s = 0.25, ^ = oo, = oo 

(b) = 0.6, e = 1, 5 = 1.25, ^ = 0.78496, ^2 = °o 

(c) k^ = 0.5, e = 0.5, s = 0.615920, ^ = 1.1315, qz = 3 

(d) k^ = 0.3, e = 0, s = 0.901890, ^ = 0.78606, ^2 = 1-1805 

(e) k^ = 0.095385, e = 0, s = 10, ^ = 0.13684, ^2 = 0.03267 




0.2 



C 









0.02 
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Fig. 8. Unbounded orbits for 

(a) = 0.99, e = 1, s = 0.2499968435, ^ = 4.6378 

n 

(b) = 1, e = 1, s = 0.25, ^ = 00 

n 

(c) k^ = 0.9998000799, e = 1, s = 0.2501 , ^ = 5.0816 

n 



Fig. 9. Trajectories of light for 

(a) f/i = 0.42922, ^(p = n/2 

(b) f/i = 0.56808, A(p = n 

(c) f/i = 0.62334, Acp = 3n/2 

(d) i/i = 0.64720, Acp = 2n 



Fig. 10. Trajectories of light that terminated 

(a) f/i = 0.8333, = 3.8345 

(b) ^ 1, (p^ = n 

(c) f/i = 10, = 0.78133 

a 




006 
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